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Abstract 

Arrow diagram (AD) method (L. Kantorovich and B. Zapol, J. Chem. Phys. 96, 8420 (1992); ibid, 
8427) provides a convenient means of systematic calculation of arbitrary matrix elements, (>I'|0|^'), of 
symmetrical operators, O, in quantum chemistry when the total system wavefunction ^ is represented as 
an antisymmetrised product of overlapping many-electron group functions, $a, corresponding to each part 
(group) A of the system: ^' = ^Y\a extended (e.g. infinite) systems the calculation is somewhat 

difficult, however, as mean values of the operators require that each term of the diagram expansion is to 
be divided by the normalisation integral S = {*| which is given by an AD expansion as well. A linked 
AD theorem suggested previously (L. Kantorovich, Int. J. Quant. Chem. 76, 511 (2000)) to deal with 
this problem is reexamined in this paper using a simple Hartree-Fock problem of a one-dimensional ring of 
infinite size which is found to be analytically solvable. We find that corrections to the linked AD theorem 
are necessary in a general case of a finite overlap between different electronic groups. A general method of 
constructing these corrections in a form of a power series expansion with respect to overlap is suggested. It 
is illustrated on the ring model system. 



1 Introduction 

For a wide range of quantum-mechanical systems, the entire electronic system can conveniently be split into 
a set of electron groups (EG) such as core and valence electrons in molecules or crystals, electrons on atoms 
or ions in atomic or ionic solids, core and bond electrons in strongly covalent materials, etc. [TJEj;2j. Similar 
partitioning ideas can also be applied to separate electrons in a cluster and environment regions to derive a 
particular embedding potential for the quantum cluster pi Bl E] . Provided that the partition scheme appHed 
to the given system is physically (or chemically) appropriate, one can assume that electrons fixed to the given 
group spend most of their time at this group. Therefore, to a good approximation, the wavefunction of the whole 
system consisting of electronic groups / can be represented as an antisymmetrised product of wavefunctions 
^i{Xi) of every individual group: 

M 

■^{Xi,---,XM) = AY[<i>j{Xj) (1) 

1=1 

where M is the total number of EG's. Here the I-th group is associated with Nj electrons which coordinates 
are collected into a set Xi = {xi, . . .xn,). Note that the single electron coordinates x = (r, a) include the spin 
coordinate a as well. The antisymmetrisation operator is defined via 



where the sum runs over all N\ permutations P of the complete group of permutations Sn, the total number of 
electrons in the whole system being N = Ni + N2 + ■ . ■ + Nm and the factor ep = ±1 according to the parity of 
the permutation P. It is assumed here that the group functions are already individually antisymmetric. 

By applying expression ^ for the operator A in Eq. (QJ , one obtains an expansion of the wavefunction via 
various functions-products. When used in calculating matrix elements (^Pl O of symmetrical (with respect to 
permutations of electronic coordinates) operators O = 0(xi, . . . ^xn) (e.g. an external field or electron-electron 
interaction), a corresponding expansion of matrix elements between various functions-products is obtained. 
Earlier attempts (see, e.g. |H1 El UHL IllL I12j l to simplify the expansion and associate each term with 
a transparent diagram have been generalised in Refs. \1?>\ I14j where the Arrow Diagram (AD) theory was 
developed. Firstly, by exploiting a double-coset decomposition of the group Sn , many terms in the expansion 
were found to be identical and the corresponding rules to number all distinct terms were worked out in a general 
case. Secondly, each distinct term in the expansion was given by a well-defined picture (arrow diagram) which 
allowed to construct simple rules to associate an analytical expression with any of the AD's. The theory was 
formulated in a very general way for arbitrary number of electron groups and explicit AD expansions were 
constructed for reduced density matrices (RDM) ^ [2| of orders one and two (RDM-1 and RDM-2). These 
enable one to calculate matrix elements of arbitrary symmetrical one- and two-particle operators and are thus 
sufficient in most cases relevant to quantum chemistry of molecules and solids (generalisation to higher order 
RDM's is cumbersome but straightforward). 

In general, the success of this technique depends crucially on the value of the overlap between different group 
functions <i>/(X/) and thus their localisation in certain regions of space (also known as structure elements 
This is because greater localisation of the group functions results in better convergence of the AD expansion, 
i.e. smaller number of terms in the expansion are to be retained. That is why the choice of electronic groups 
and appropriate regions of their locaHsation is so important for this method to work. 

Although the AD theory formulated in \1?>\ ITl] is applicable to any system consisting of arbitrary number of 
overlapping electron groups, one can not directly use it to describe an extended (in particular, infinite) system. 
Indeed, to calculate an observable 

(vl/|vl/) 

of the operator O, in the case of an extended system one has to calculate the ratio of two AD expansions, one 
arising from the matrix element (4'|0|4') and another - from the normalisation integral S = {^\ ^) . Both 
expansions contain very many terms; in fact, the number of terms is infinite in the case of an infinite system 
(e.g. a solid). 

It was argued in pj] that the normaHsation integral S tends to infinity with the number of electronic groups 
M . This assertion made it possible to prove the so-called linked AD theorem QJ]. It states that the observable 
O — (^'1 O l^*)^, i.e. it is equal to the AD expansion for the matrix element (^'|0|\I') of the operator O, 
in which only linked (connected) ADs are retained (subscript "c"), and the normalisation integral should be 
dropped altogether. This expression is very attractive as it allows one to consider a single AD expansion for 
arbitrary matrix elements of operators. 

It is the main objective of this paper to show that the Hnked AD theorem is correct only approximately. 
When the overlap between different group functions becomes larger, corrections are to be introduced. We show 
this by first considering in detail a simple "toy" model of a ID ring for which the exact analytical solution is 
possible. This allows us to investigate in detail the limiting behaviour of the AD expansion for this system when 
the number of groups in the ring tends to infinity. Then, we suggest a general method of calculating the mean 
values, Eq. lO, in a form of corrections to the Hnked AD theorem. 

The paper is organised as follows. In the subsequent section we shall briefiy review the AD theory and 
introduce the linked AD theorem for the normalisation integral and RDM-1. The consideration of higher order 
RDM's is analogous. In section 3, we consider in detail the toy model and obtain an exact expression for RDM-1 
in the Hmit of an infinite ring size. This exact results will then be compared with that obtained using the Hnked 
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AD theorem. A systematic way of constructing corrections to the Hnked AD theorem is then suggested in 
section 4. 



2 The Arrow Diagrams theory 

First we consider a symmetric group Sni ■ Its elements permute electronic coordinates belonging only to the I-th 
group. Joining together all such groups we obtain a subgroup Sq = Sni U U • • • U Snm of Sn- Any element 
of 5*0 only interchanges electronic coordinates within the groups, i.e. performs only intra-group permutations. 
The algebraic foundation of the AD theory is based on the double coset (DC) decomposition of the complete 
symmetric group Sn with respect to 5*0, which enables to single out all inequivalent intergroup permutations 

^ TtYI fJ-qTSoPqTSo (4) 
iVn — 
" qT 

Here the sum runs over all distinct types, q, of operations for intergroup permutations as well as all distinct 

ways, T, of labelling actual groups (EC's) involved in it, and Nq — Ni\N2l A^m! is a numerical factor. 

PqT is a DC generator, which, involves only inter-group permutations and, in general, can be constructed as a 
product of some primitive cycles each involving no more than one electron from each group. 

Each such a cycle is represented as a directed closed loop connecting all groups involved in it. Thus, 
in general, each permutation PqT can be drawn as an arrow diagram containing a collection of closed directed 
loops. These loops may pass through the given group several times (depending on the number of electrons of the 
group which are involved in the permutation Pqx)- If the groups involved in the diagram cannot be separated 
without destroying the directed loops, the AD is called linked or connected. If, however, this separation is 
possible than the diagram is called non-linked or disconnected and it can be represented as a collection of linked 
parts. The decomposition coefficients fiqT can be calculated merely by counting arrows entering and leaving 
each group in the AD. 

By writing the antisymmetriser A of Eq. ^ using the DC generators, it is possible to obtain the diagram 
expansion of the normalisation integral jl4] : 

s = (a<^ 



A^)=AY^ eqfiqT PqT $) (5) 

where A = Nq/N\ and $ — Ylj is a product of all group functions. Eq. Q gives an expansion of S in 

terms of diagrams identical to those for the permutation group Sn of Eq. Q . The matrix elements ( $ 



qT 



are represented as a product of RDM's of electronic groups involved in the permutation PqT integrated over 
corresponding electronic coordinates. If each of the group functions <&/ is given as a sum of Slater determinants 

'%T 



with molecular orbitals expanded via some atomic orbitals (AO) basis set, then the matrix elements PqT ^ 
are eventually expressed as a sum of products of simple overlap integrals between the AO's of different electronic 
groups (see, .e.g. Ref. ^|). It is important to note here that a total contribution of a non-Hnked AD is exactly 
equal to the product of contributions associated with each of its linked parts |14j . 

For the following, it is convenient to sort out all the terms in the diagrammatic expansion by the number of 
groups involved in the AD's. Getting rid of the common factor A, we can write: 

^=^ = 1+E E SK{A^,A2,■■■,AK) (6) 

K=2 Ax< - <Ak 
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where SxiAi, A2, ■ • • ,Ak) is the contribution of all AD's that correspond to intergroup permutations among 
K different groups with the particular labelUng ^1,^2,... ,Ak. The unity in equation above corresponds to 
the trivial permutation. 

For convenience, the expansion in Eq. ^ will be referred to as the normalisation-integral expansion. In 
principle, in this expansion all groups of the entire system participate. It is also found convenient to introduce 
a derivative object, S{T), which is obtained from the above expansion by retaining all the AD's which are 
associated only with the groups from a finite manifold T — {Ai, . . . ,Ak} (i.e. containing groups with labels 
Ai, . . . , Ak). If T comprises the whole system, we arrive at S of Eq. (j^J. We shall also introduce a manifold 
[T], i.e. an artificial system which is obtained by "removing" from the entire system all groups comprising the 
set T. Then, S'([T]) (we shall also use a simpler notation S[r] when convenient) is obtained from Eq. Q by 
retaining only AD's in which any of the groups belonging to T is not present; in other words, ^([T]) is obtained 
from S by assuming that in any of the AD's all the overlap integrals involving groups from the manifold T are 
equal to zero. Obviously, S can be considered as a particular case of S[r] when the manifold T is empty. We 
shall occasionally also use the notation S[Ai, . . . , Ak] for S\T], where T = {Ai, . . . , Ak}^ if we want to indicate 
explicitly which particular groups are excluded. 

The unnormalised RDM-1 of the whole system. 



can also be written as a matrix element of a certain symmetrical one-particle operator 0|2lE]. Therefore, as 
in the case of the normalisation integral, by inserting the AD expansion for the operator A, one obtains the 
corresponding AD expansion for the RDM-1 pi]. It can be constructed by considering the AD expansion of S 
and then modifying each diagram by placing a small open circle, representing the variables {x,x'), in either of 
the three following ways: (i) on a group not involved in the diagram; (ii) on a group involved in the diagram, 
and, finally, (iii) on an arrow. Thus, each AD in the S expansion serves as a reference in building up the AD 
expansion for the RDM-1. Since the contribution of any diagram, consisting of non-linked parts, is equal to 
the product of contributions corresponding to each of the parts, it can be shown ^2] that, in general, the AD 
expansion of the RDM-1 can be represented as 



where p^(Ai, . . . , Ak \\ x; x') is the sum of contributions of all AD's with an open circle which are constructed 
using the particular group labelling ^1, ^2, • ' ' i ^k- The sum over t takes account of the two positions of the 
open circle on the reference AD: on an arrow and on a group. Note that each AD with the open circle is 
necessarily a linked (connected) AD. One can see from the above equation that each AD with an open circle, 
constructed from groups T — {Ai, . . . , Ak}, is multiplied by the sum of all possible normalisation- integral AD's 
(i.e. AD's without the circle), S[Ai, . . . , Ak] = S[T], constructed using the rest of the system [T]. 

The RDM-1 introduced by Eq. ijjjl is not normalised to the total number of electrons in the system since the 
wavefunction 5" constructed using group functions, $/, is, in general, not normalised to unity: S — {'^\ ^) ^ 1. 
Therefore, the true RDM-1 should be calculated according to Eq. iPJ as |17) : 




(7) 




M 



K=l Ai<...<Ak 




M 



(9) 



where 



f[T] - f[Ai,...,AK] 



S[Ai,.^.,Ak] 
S 



(10) 



4 



(a) 



(b) 



Figure 1: The ID 'toy' model: (a) a HF ring consisting of equidistant one-electron orbitals ipi, ip2, ■ • •, and 
(6) the corresponding chain model in which there is no overlap between the first and the last orbitals in the 
ring. 

are numerical pre-factors. These depend only on the chosen set of groups T. Since the normalisation integral 
S is represented via an AD expansion, Eq. J^l, one can see that the pre-factors are given as a ratio of two 
AD expansions, each containing a very large (infinite) number of terms for a large (infinite) system. Thus, it 
follows from the last passage in Eq. that the true RDM-1 is represented as a sum of all linked AD's with 
the open circle, p*j({Ai, . . . , Ak \\ x;x'), multipHed by numerical pre-factors, f[Ai,...,AK]- The calculation of the 
pre-factors poses the main problem in applying the AD theory to infinite or even large system (for small systems 
all AD's can be accounted for explicitly and thus p can easily be calculated, at least in principle). 

It was argued in ^2] that in the limit of an infinite system the pre-factors f^r] tend exactly to unity for any 
choice of the groups T — {Ai, . . . , Ak}- As a result, it was argued that the true RDM-1 can be represented as 
a single AD expansion containing only AD's with the open circle which, as was mentioned earner,^are all Hnked 
AD's. The main argument put forward to prove this so-called 'Hnked-AD theorem' was that any S[Ai, . . . , Ak] 
(and, in particular, S) tends to infinity when the number of electronic groups M tends to infinity. 

The main objective of this paper is to show that the situation is more subtle. Although the argument about 
the limiting behaviour of S[T] put forward in ^J] may seem quite plausible, in some cases, as will be shown in 
the next section by considering an exactly solvable 'toy' model, S[T] may have a zero limit. The main reason 
for this is that contributions of different AD's in the expansion have alternating signs, the point overlooked 
previously. As a result, the pre-factors /[t] do not necessarily tend to unity, but may take a different value 
depending on the given system. 

3 ID toy model: a Hartree-Fock ring 

Let us consider a ring of M equally spaced one-electron groups as shown in Fig. n](a). Each group is described 
by a single real s-type normalised wavefunction ipi{x) (i = 1, 2, . . . , M) localised around the group centre. All 
groups are identical, i.e. every localised function ipi{x) can be obtained by a corresponding translation of the 
neighbouring function. Note that N — M for this model and the ring wavefunction of Eq. (QJ represents 
essentially a single Slater determinant 5" oc det{V'i(xi) • • ■ ^m{xm)}- Hence, this theory should be equivalent 
to the Hartree-Fock (HF) method. 

We assume that only neighbouring functions overlap. Due to the symmetry of the model, only one parameter. 
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namely the overlap integral cr = J ipi{x)'ipi+i{x)dx between any neighbouring groups, is required to characterise 
the system; the actual spatial form of ipi{x) is not actually needed. 

We shall show that the RDM-1 of this system can be calculated analytically for any M. Most remarkably, 
it will be shown that the M ^ oo limit is also analytically attainable. Using this result, we will analyse the 
linked- AD theorem. We shall find that some corrections are required. 

Two independent methods will be used to consider the toy model. 



3.1 Finite ring: method based on the inverse of the overlap matrix 

It is first instructive to perform the calculation using a method which is completely independent of the AD 
theory. This is indeed possible since an expression for the RDM-1 of any HF system can be written explicitly 
via orbitals and their overlap matrix Q]. Taking into account that only overlap between nearest neighbours 
exists, we obtain for the ring electron density (the diagonal elements of the RDM-1): 

N 



N 

E 

1=1 



Here is the inverse of the overlap matrix 



(11) 



1 

(7 



(7 
1 



1 



(12) 



NxN 



between the ring orbitals. Note that all the omitted elements in the above expression are zeros. This convention 
is adopted hereafter. Then, the normalisation integral: 

G^^"^ = detS (13) 

is simply given by the determinant of the matrix S of Eq. II12II . To avoid a possible confusion with the overlap 
matrix, we use for the normalisation integral a different notation here, G^*"^. 

For the convenience of the following calculation, we also introduce a similar chain system shown in Fig. 
(6) in which the overlap between the groups 1 and M is equal to zero. The normalisation integral for the chain 
is 



G Chain 
N 



cr 
1 



1 
a 



(14) 



NxN 
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Owing to the symmetry of the above introduced determinants, it is possible to derive the following recurrence 
relations: 

G^"9 = G%'IT - 2a^G%'lT + 2 (-1)^+' (15) 

G Chain r^Chain 2r^Chain (^c^\ 

The first one is obtained by opening the determinant of the matrix Ijl2ll along its first row (or column) , while 
the second relation is obtained in the same way from Eq. ifnji . Combined with the obvious 'initial' expressions 
^Ring ^ QChmn ^ ^ QChatn ^ 1 _ cT^, the above relations allow one to calculate the GC/mm QRing 
any value of N . For instance, we obtain, 

^mng ^ 1 _ 6(72 + - (17) 
G'i^""' = 1 - 5(7^ + Qa'^ - (19) 

QChain ^i_Q^2 ^ -^Q^4 _ 4^6 (20) 

Note the alternating sign in every next term in the expressions above. Also, in agreement with one can see 
that numerical pre-factors in each term grow with increase of N. 

It is now straightforward to investigate the limiting behaviour of these quantities when the number of groups 
tends to infinity (the N ^ oo limit). We find, using a simple numerical calculation, that both of them tend 
to zero with increase of N for any value of cr < 0.5 (it will be clear later on that the limit does not exist for 
a > 0.5). Since both G^'"'™ and G^J"'^ correspond to the normalisation integrals for the two systems, one 
can see that this limiting behaviour is quite different from what should have been expected from the analysis 
performed in Ref. This is explained by the fact that both limAf^+oo G^*"^ and limAf^oo QCham ^^^^ given 

by infinite series with alternating terms. The limiting behaviour of each of the series is not straightforward in 
spite of the fact that the pre-factors to each term cr"(n = 2, 3, . . .) tend to infinity when N ^ oo. 

In order to calculate the RDM-1 lfTT|l . one also need the corresponding elements of the inverse of the overlap 
matrix. Using the explicit structure of the matrix S given in Eq. lfT2ll . one obtains: 

^ ^'-^ detS G^'"^ 
_ _ _ -crG^'^J'^'^ + (-1)^+^ cr^~^ 

" \^ ~ dg^S ^ QRing y'^'^l 

where Gi^i and C^.i+i are cofactors of S^^i and ^i+i.i respectively. 
Using Eqs. Ij21ll and (I22II . we obtain for the electron density ljll|l : 

p{x) = 




N o // \N\ N 



/ i=l \ / i=l 



The above expression allows one to calculate the electron density of the ring with arbitrary number of groups. 
The N —> oo limit will be considered separately in section 13.31 
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Figure 2: The AD expansion for the ring normalisation integral, S. The ring consists of a finite number N of 
one-electron groups. 




Figure 3: The AD expansion for the diagonal elements of the ring RDM-1, p{x) = p{x;x), Eq. 
3.2 Finite ring: method based on the AD theory 

We shall show in this section that the same expressions as derived in the previous subsection can also be obtained 
using a much simpler algebra of the AD theory. ^ 

To this end, we first consider the normaHsation integral S = G^™^ of the ring. Since all groups contain 
a single electron and only nearest groups have non-zero overlap, the AD expansion contains essentially only 
bubble-like AD's between two adjacent groups (including all non-linked AD's constructed out of them) and two 
AD's connecting all groups (a TV-vertex polygon) with opposite direction of arrows, as depicted in Fig. |21 An 
AD expansion for the G^'"*™ looks similarly. There are two differences: (i) bubble AD's containing the groups 
1 and N are missing, and (ii) it does not contain the two polygon diagrams. Recall that there is no "connection" 
between the groups 1 and N in the chain. 

An AD expansion for the electron density can also be written explicitly as shown in Fig. El It contains three 
terms. The first term is associated with an open circle AD, associated with, say, group i, multiplied by all bubble 
AD's (shown in the square brackets) which are constructed out of the rest N —\ groups — l,i-fl,...,iV. 

The contribution of the open circle pf\x) = tpfix), while the bubble AD's in the square brackets all amount to 
QChain ^g^jjj g^j,g ^YiQ same for each ?), since they represent all AD's for a chain with A^— 1 groups. Finally, one has 

to sum over all values of? (all positions of the open circle). Hence, the first term contributes G^'^f"' X^iLi Pt'\-'^) 
to the density p(x). 

The second term in Fig. Elcontains two bubble diagrams with an open circle on both of its arrows, multipHed 
by all possible bubble AD's constructed out of the rest N — 2 groups. Assuming that the bubble AD is 
taken between the groups i and i + 1, we obtain the contribution pf^^'f'^^x) = —a')pi{x)')pi+i{x) for it. Since 
the contribution of all groups — l,i + 2, . . . , N in the square brackets in the second term is equal to 

QChain (^hich is the same for any pair i, i + 1), the final contribution of the second term in Fig. becomes 

SG^'^fj" J2iLi pfi+i'^\^)^ where we have also summed over all pairs, i.e. all possible positions of the open circle 
on the bubble AD's (recall, that the bubble AD's can be taken only between nearest neighbours in our model). 
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Finally, the last two terms in Fig. Elare represented by polygons with an open circle. Each such a diagram 
contributes p^f°l\' = {—a)'^~^ 'ipi{x)ilji+i{x) if the open circle is positioned on the arrow connecting groups i 
and i + The total contribution of the polygons is obtained by summing over all possible positions of the circle 
and multiplying by a factor of two since each AD accepts two directions of the arrows. 

Summing all three contributions, we arrive at the following expression for the unnormalised electron density: 

N N N 

m = G%'iT E (^) + E f^^^ (^) + 2 E (^) (24) 

i=l i—1 2—1 

Dividing this expression by G^*"^, one obtains the normalised electron density p{x), which, as can easily be 
seen, appears to be identical to that of Eq. obtained in the previous section using Slater determinants. 

To finish the proof of the equivalence of the two methods, we only need to re-derive the recurrence relations 
lfT5|l and ifTfiji using exclusively the AD theory. 

Consider first the expansion for G^'"'™ which contains only zero, one, two, etc. bubble diagrams (see the 
bubble-containing terms in Fig. I^J. Recall that the chain is broken between groups 1 and N. Let us "fix" group 
1. Then, G^'"'™ can be written as the sum of two contributions: (i) due to all diagrams involving group 1 
and (ii) due to the rest of them which do not involve it. The first contribution is simply a single bubble AD 
between groups 1 and 2, equal to — cr^, times all the bubble AD's made of all the other groups, i.e. the whole 
contribution is (— cr^) G^'^""- The contribution of all AD's which do not involve group 1 is simply equal to 
can immediately recognise the recurrence relation 

Eq. ltT5|l is proven in a similar way. Again, we first "fix" group 1. Then, the sum of all AD's shown in Fig. 
His equal to the sum of three contributions: (i) G^'^'Y" due to all AD's which do not involve group 1; (ii) all 
the bubble AD's involving it, given by 2 (-cr2) G^'L'^a" analogously to the case of the chain discussed above (the 
factor of two arises due to the fact that one can construct two bubble AD's with group 1, namely those involving 
group pairs 1, and 1,2); (iii) the contribution —2 (— a)^of the two polygon AD's. Summing all three terms, 
one arrives exactly at Eq. Ijl5|l . as required. 



3.3 Infinite ring 

Our task now is to calculate the pre-factors (in the round brackets) in the expression for the density of Eq. 
H23|l . To do this, it was found convenient to introduce the following quantity: 



Pn=-'\t^\ (25) 



Using explicit expressions for the Gj''"*™ and Gf"^'", we get P2 = cr^/ (l — f^) • Using the recurrence relation 
for G^''"™, a very simple recurrence relation for P/v can also be derived: 

= , p (26) 



The introduced quantity is very useful in two respects. Firstly, the first two pre-factors in Eq. I|23|l can be 
directly expressed via it: 

^^-1 ^ t f27l 

QlUng \-2PN-l+2aXN-l ^ ' 
r^Chairi ^ — 2 r> 

a Pa._i ^28) 



^mng I _ 2P^_^ + 2axN-i 
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where ^ 

Secondly, it has a finite N ^ oo limit. Indeed, assuming that such a limit exists, we obtain from Eq. Ij26|l 
a simple quadratic equation Poo(l — ^oo) = tr^, where Poo — ^t^^n^qo Pn- The roots of this equation are 
^ (l ± Vl — 4(t2) . To choose the correct sign, we consider the case of cr = 0, when G^'"^™ = 1 for any N, and 
thus Pn = 0. Only the root with the minus sign, 



1 - Vl - 4cr2 

^oo ^ (30) 

satisfies this condition. It also appears that a definite real limit exists only for <J < \ which corresponds to a 
not very large overlap. This restriction is a consequence of the nearest-neighbour approximation adopted in our 
toy model. Note in passing that Pn > for any TV. 

Finally, to finish the calculation, we have to consider the limiting behaviour of the ratios xn and 



which enter Eqs. Ij27|) . (I28II and (I23|l . respectively. To calculate the limj\j^oc Xn, we use the definition Il25|l of 
Pat to write: 

Chain V / /^Chain I ^ \ /^Chain 



• • • ^3 V ^oo / 

since Pn < Poo for any finite value of N and any < cr < ^. Since G^'"^™ = 1 — > and any of the 
quantities Pk are positive, one can write 



< --TTT— < 



r^Chain /"'Chain 



Po_ 

a 



Because < 1 for any < cr < i, we obtain limjv^oo ("^)^ = and then 



Hm Z,. ■ I = (32) 



Hence, xn of Eq. Ij29ll has a zero limit and disappears in the ratios II27II and Ij28|l . Consequently, ^at ^ too 
(see Eq. (jHH), so that the third term in the right hand side of Eq. lf23jl also does not contribute when the limit 
of an infinite ring is taken. 

Finally, combining all contributions, we obtain the following exact result for the infinite ring: 



N N 

p{x) = ^ f^f^ix) + 2^ /[,,+i]p|;:tf ^^(-) (33) 



i=l i=l 

where the corresponding expressions for the two pre-factors are: 



^[''^+^1 = i (tt^ " ' ^ ^^^^ 
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Figure 4: Comparison of the exact expression for the pre-factor /[i]((t), Eq. with its series expansion, Eq. 
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Figure 5: Comparison of the exact expression for the pre-factor fii.2]{<^), Eq. . with its series expansion, 
Eq. ggi. 
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The above obtained expression for the normahsed density can be compared directly with the general formula 
((HJ. According to the linked- AD theorem both pre-factors are equal to unity for any a. We see, however, 
that is not the case in general. The actual pre-factors /[ij and are plotted in Figs. 21 and El respectively, 

for the whole range of the overlap integral a by solid Hues. 

One can clearly see that for < ct ^ 0.3, i.e. for the small and intermediate overlap, both pre-factors 
are approximately equal to unity, which is the value predicted by the linked- AD theorem. However, when the 
overlap between neighbouring group functions is larger, 0.3 ^ a < 0.5, both pre-factors shoot up to infinity. 
Note that the singularity at cr = 0.5 is somewhat artificial and is removed if the nearest-neighbour approximation 
for the overlap is lifted. 

Thus, we conclude, that the linked- AD theorem works well in the region of not very large overlap, i.e. when 
the group functions are sufficiently localised. If the overlap is large, then corrections to the linked- AD theorem 
are necessary. A general method for building up such corrections will be suggested in the next section. 

4 Corrections to the linked- AD theorem 

The idea of the method to be proposed in this section is based on the power series expansion for the pre-factors. 
In order to illustrate the general method, it is instructive to consider first a simple example of the familiar ID 
toy model. 

4.1 Series expansion for the ID toy model 

Let us consider the ID toy model (an infinite ring) for which the exact solution is known. Fix group 1. Then, 
the normalisation integral for the whole ring, S, can be constructed as a sum of three terms: (i) all AD's of the 
rest of the system, S'[l], where T — [1] corresponds to the ring with group 1 "removed", i.e. this group does not 
have any overlap with its neighbours; (ii) a bubble AD between groups 1 and 2 (contributing —cr^) multiplied 
by all possible bubble AD's due to all other groups, S[l, 2] (as usual, [1, 2] denotes the system in which groups 
1 and 2 are removed), and similarly (iii) a bubble AD between group 1 and the last group in the ring (the 
other neighbour of group 1) times all the bubble AD's of the rest of the system. The last two contributions 
are identical. Note that we do not consider the polygon here as we know f section l3.3|l that in the limit of an 
infinite system its contribution vanishes. Thus, one can write: S = S[l] - 2a'^S[l,2]. Dividing both sides of 
this equation by S, we obtain an equation connecting pre-factors: 

1 = /[I] - 2(7V[i,2] (36) 

Repeating the above procedure for the system [1] and fixing group 2, we similarly obtain: S[l] — S[l,2] — 
a'^S[l, 2, 3], where [1, 2, 3] denotes a system in which groups 1, 2 and 3 are removed. In fact, the above equation 
also follows from the recurrence relation Ijltill . Dividing the last equation by S, we obtain another relationship 
between pre-factors: 

/[I] = /[1, 2] - 0-^/[l,2,3] (37) 

Continuing this procedure, one can write an infinite series of relationships in which every time a new pre-factor 
appears, for instance: 

/[1, 2] = /[1, 2, 3] - 0-^/[l,2,3,4] (38) 
/[1, 2,3] = /[1, 2, 3,4] - O-^ /[1, 2, 3, 4,5] (39) 

Next, we assume that the pre-factors are well-defined functions of the overlap, a. Hence, they all can be 
expanded in a power series: 

oo 

/[T]=l + E«iri'^" (40) 

)i=i 
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where T is either 1 or 1, 2 or 1, 2, 3 and so on. Note that the expansion starts from the unity for any [T] since 
fyr] = 1 for zero overlap ((7 = 0). Substituting these expansions into Eqs. Ij36|l - (j39ll . one can recursively recover 
the expansion coefficients up to a certain order by comparing terms of the same power of a. Indeed, it follows 
from Eq. ^ that af' = 0, a|,^' = 2 and aL^' = 2a[,^l^] for any ti > 3. From Eq. ^ we get: = a^^^' = 0, 

a2^'^' = 02^' + 1 = 3 and a[i' — ftn'^' + ~ ^ for any n > 3. Similarly, we obtain from Eq. Ij38|l that: 

[1,2,3] [1,2] [1.2,3] [1,2] , -, . , [1,2] [1,2,3] , [1,2,3,4] „ r o n u- • j-u u 

a!^ = a!^ = (J, 02 = 02 + 1 = 4 and an ' — an ' + ^^-2 = U tor n> 6. Combining the above 
recursive relationships, new coefficients can be obtained, e.g. flg^'^' = 0, 04 = 10, 03^' = 0, O4 ' = 6, Og^' = 0, 
aj;^' = 20. Using Eq. Ij39|l enables calculation of more coefficients. This way we obtain the first several terms in 
the expansion of the two pre-factors that are required for the ID toy model: 

/[I] = 1 + 20-2 + Qa'^ + 20(T^ + 70(T^ + . . . (41) 



/ji 2] = 1 + 3(7^ + lOfj" + 35cr^ + . . . (42) 

It can easily be seen by expanding in a power series the exact expressions l(3ijl and l(35)l for the pre-factors, that 
the above expansions are indeed correct. To obtain more terms in the expansion by using this method, one has 
to consider higher order pre-factors. The expansions for /[ij and /[i,2] obtained above are compared with the 
exact expressions of Eqs. Ij34|l and Ij35|l . respectively, in Figs. 01and|3 One can see that in both cases reasonable 
approximations to the exact values of the pre-factors are obtained using the power series expansion in the cases 
of small and intermediate overlap. It is also apparent that one has to go to higher orders in the expansions 
when the overlap is much larger. In particular, the power series expansions do not show sharp increase of the 
pre-factors closer to the critical overlap a = \- 

4.2 General method 

The idea of the series expansion method of the previous subsection can actually be generafised for an ar- 
bitrary system containing groups with any numbers of electrons, i.e. for the most general wavefunction of 
Eq. (QJ. Indeed, what is needed is a way to connect the pre-factor /[^-j = S\T]IS associated with some set 

T = {^1,^2, . . . , Ak} of groups (not necessarily nearest neighbours) with the pre-factor /[y] = S\T']/S asso- 
ciated with a smaller set \T'] = [T IJ AT] = \T + AT] . The set T' contains L additional groups forming a set 
AT = {Bi, i32, • ■ • , -Sl}, see an illustration in Fig. El 

The corresponding relationship between the two pre-factors can be^obtained as follows. Consider the sum 
of all possible AD's which make up the entire normalisation integral S\T] for the subsystem [T]. From all of 
them we separate out those which are related to the set [T -I- AT] : 

S'iT] = S'(AT)S'[T-|- AT] + ^ ^ i5(TiT2)5(AT - Ti)S'[T AT T2] (43) 

TiCAT T2<Z[T+AT] 

The first term in the right hand side contains all the AD's that can be constructed out of the groups of the set 
AT, given by S'(AT). These should be multiplied by all AD's constructed out of the rest of the groups, i.e. by 
all groups from the manifold [T -I- AT]. The second term contains all the AD's which are formed by the groups 
belonging to both sets AT and [T-|- AT]. Indeed, the double sum picks up subsets Ti and T2 from the manifolds 
of groups AT and [T -I- AT], respectively. The sum of all possible AD's constructed using every group from Ti 
and T2 is denoted ^'(TiTj). Note that AD's in D{TiT2) may contain non-linked diagrams as well; however, in 
this case each of non-linked AD's must contain all groups from Ti and T2. Every AD in D{TiT2) is multipHed 
by all possible AD's formed out of the rest of groups of the two sets: S{AT — Ti) corresponds to the AD's 
made of all groups in the set AT — Ti (created by removing groups belonging to Ti from the set AT) , whereas 
S[T + AT + T2] is obtained by all AD's left after "removing" the set T2 from the set [T + AT]. Eq. ^ is 
schematically illustrated in Fig. d 
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Figure 6: A schematic of an infinite 2D lattice of electronic groups. The solid crosses indicate the "removed" 
groups that comprise the set T = {a, 6, c, d, e}, whereas dashed crosses indicate additional groups to be "removed" 
which form a set AT = {/, g. ft,}. 
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Figure 8: AD's representation oi D{AB) of Eq. lUfill . 

Dividing each term in Eq. by S, we obtain a recurrence relation for the pre-factors sought for: 

/[T] = 5(AT)/[T+AT] + E DiTiT2)S{AT~T,)f[T+AT+T,] (44) 

TiCAT T2C[T+AT] 

It is seen that /[t] is represented as a Hnear combination of pre-factors corresponding to smaller systems: if 
Ht is the number of groups in the manifold T, then we can write that ^[t] > ^[t+at] > ^[t+at+Ts] ■ Eq. ijiiji 
is a generalisation of any of the Eqs. - written above for the toy model. 

Equations like the one written above allow one to obtain the necessary series expansion in the general case. 
To this end, we attach to every diagram a factor cr" with the power n being the order of the AD, i.e. the number 
of arrows it contains. The expansion we obtain will be with respect to a that will be set to unity at the end of 
the calculation. Note that this method does not take into account an additional factor, which is relevant in the 
actual calculation, that the contributions of the AD's also depend on the distances between different groups. In 
the method we propose an expansion with respect to overlap integrals will be obtained. 

To illustrate how the method works, we shall obtain a few first terms in the expansion for /[^j , /[a,s] and 
flA.B,c]- Noting that S{A) = 1, we start, similarly to section ITTl with T = (empty) and AT = A in Eq. ijiill : 

1 = f[A] + D{AB)f[AM + J2 D{ABC)f[A,B.c] + J2 D{ABCD)fyA,B,c,D] + ■■■ (45) 

Be[A\ {B,C}£[A] {B,C,D}e[A\ 

where the curly brackets under the summation sign indicate that we sum over all sets of the groups containing 
two, three, etc, groups in them irrespective of their order in the sets; further, 

b{AB) = D2{A, B)a^ + Di{A, B)a^ + ... (46) 

D{ABC) = D3{A, B, C)a^ + D^iA, B, C)a^ + D^iA, B, C)a'' + ... (47) 

contain all AD's composed of groups {A, B} and {A,B,C}, respectively. For convenience, contributions of 
different orders with respect to a have been presented separately by the quantities like Z3„(T). In particular, 
D{AB) contains a sum of bubble AD's which only have even number of arrows as shown in Fig. |H| The 
expansion is finite: it runs until the number of electrons in either of the groups is completely used up. For 
instance, i^a group A contains only two electrons, there will only be the first two terms in the expansion. 
Similarly, D{ABC) contains a finite sum of all AD's made out of groups {A,B,C} as shown in Fig. El (note 
that both directions of arrows, where appropriate, are assumed in this expansion). The expansion starts from 
the AD containing 3 arrows, the next term contains 4 arrows, then follows the term with 5 arrows, etc. For 
small values of n there might be only few AD's contained in the given £)„(r); however, the number of the AD's 
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Figure 9: AD's representation of D{ABC) of Eq. ll47jl . 

may grow quickly as the order n becomes larger. The only exception is D{AB) in which case in any order n 
(that is always even) there is exactly one AD in Dn{T). 

Following the same logic, consider Eq. ijiiji for T = A and AT — B. This gives: 

I[A] = hA,B]+ E D{BC)JyA,B,c]+ D{BCD)f[A,B,c,D] + --- (48) 

Ce[A,B] {C,D}e[A,B] 

If r = {A, B} and AT = C, then we also have: 

f[A,B] = flA,B,C] + ,B,C,D] + ^ D{CDE)f^A 

,B,C,D,E] + ■ • ■ (49) 

De[A,B,C] {D,E}e[A,B,C] 

This process can be continued: at every next step one chooses the set T + AT of the previous step as the set 
T, and one group from the new [T] is chosen as the set AT. 

Expanding all the pre-factors in terms of a (see Eq. ijifijl l and comparing terms with the same powers of a, 
we obtain after setting a = 1: 

/[Al = l- E D2{A,B)- MAB,C)+ Y iD2{A,B)]' 

Be[A] {B,c}e[A] Be[A] 

+2 Y D2{A,B)[D2{A,C) + D2{B,C)]-YD'ii'^'By ^ D4{A,B,C)- ^ DM,B,C,D) 

{B,C}e[A] Be[A] {Bfi}e[A] {B,C,D}g[A] 

+ {Dz{A.B,C)[AD2{A,B)+D2{B,C)]-D^{A,B,C)}- Y D^iA, B,C, D, E) 

{B,C}e[A] {B,C,D,E}e[A] 

+ Y {(iD3iA,B,C)[D2{A,D)+D2iB,D)]+3D3{B,C,D)D2{A,D)-D5iA,B,C,D)} + ... (50) 

{B,C,D}£[A] 

The expression for the pre-factor /[^j was obtained up to the 5-th order with respect to overlap and with the 
terms ordered appropriately. Expansions for and /[a,b.c] a-re obtained up to the fourth order with respect 

to overlap using the same equations l(l5ll . and l(i9ll . Because in the final expression for the RDM-1 the 
pre-factors are multiplied by the Hnked AD's (with the open circle) which are also of some order with respect 
to the overlap, third and second order expansions are needed for fiA,B] and /[a,b.c]) respectively, if the same 
order expression for the RDM-1 is to be derived. Then, one obtains: 

/[AS] = 1-^2 (AS)- Y [D2{A,C)+D2{B,C)]- Y DM,B,C) 

Ce[A,B] C£[A,B] 

- Y [D^{A,C,D)+D^{B,C,D)] + ... (51) 

{C,D}e[A,B\ 
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Figure 10: The first three terms in the AD expansion for the RDM-1. 

f[A,B,c] = l-[D2{A,B)+D2iB,C)+D2iA,C)]- [D2{A, D) + D2{B , D) + D2{C\ D)] + . . . (52) 

DelA,B,C] 
IT] 

It is seen that in all orders the first-order term vanishes, i.e. a[ — 0. 

Note that the choice of manifolds T and AT is somewhat arbitrary. Any choice results in the same expansions. 
However, the method outHned above is, in our view, the simplest one which leads to a systematic derivation of 
the expansion coefficients for the pre-factors of increasingly higher orders. 

5 Conclusions 

In this paper we have re-examined the "Hnked-AD theorem" of Ref. ^2] formulated for extended (infinite) 
systems containing a very large number of electronic groups. For example, in the case of the electron density, 
the theorem states that the pre-factors, f^T], to the linked AD's with an open circle (involving groups of the 
set T) are all equal to unity. By analysing in detail an exactly solvable system (ID Hartree-Fock ring with a 
single electron in each group and the nearest neighbours non-zero overlap), we find, however, that this theorem 
is only approximately valid when the group wavefunctions are strongly localised. When this is not the case, so 
that the overlap between neighbouring groups is significant, corrections to the theorem have to be applied and 
the pre-factors may strongly deviate from the value of unity. 

To obtain such corrections, a general method based on the expansion of the pre-factors in power series with 
respect to overlap between different groups, has been suggested. To illustrate the method, we have obtained in 
detail the corrections for the first three pre-factors in the expression for the RDM-1, namely: (i) f[A] (associated 
with the circle AD); (ii) f[A,B] (a bubble AD with an open circle on an arrow) and (iii) f[A,B,c] (a triangle AD 
with an open circle on an arrow), as shown in Fig. Other pre-factors can be obtained along the same lines. 
Note that expressions for the pre-factors do not depend on the particular type (topology q, see section of 
the AD with an open circle; they only depend on the particular groups (the set T) involved in the diagram. 
Due to this, the same pre-factors will appear in the expansion of the RDM-2 as well. Several first terms in the 
expansion of the latter are shown in Fig. El The filled and open black circles designate the variables {x, x') 
and {y,y') of the RDM-2 p {x,y \x' ,y') , respectively (see |17l 114) for more details). 

Therefore, an energy expression up to arbitrary order with respect to overlap can finally be obtained taking 
into account the described corrections to the linked- AD theorem. It is different from that presented previously 
[Tt] only by the appearance of the pre-factors, so that we do not show it here in detail. 

It is relevant to note that an expansion of the electron density with respect to the overlap in the Hartree-Fock 
case (with no electron correlation within groups included) is well known [J. It is obtained from Eq. ifTTll by 
representing the overlap matrix S as S = 1 -I- A and then expanding = (1 -I- A)~^ with respect to the 
matrix A; the latter contains only overlap between different groups and zeros along the diagonal. The method 
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p(x,y|x,'y')= 



+ 7 f [A,B,C] 

A,B,C 




+... 



C B 



Figure 11: A few first terms in the AD expansion for the RDM-2. 

developed in this paper in section l4!2l generalises this expansion method to the case when group functions are 
linear combinations of Slater determinants, i.e. when intra-group correlation effects are accounted for in each 
group. However, it is also known (see, e.g. pHI) that the overlap expansion of the electron density in the HF 
case diverges if the overlap between groups is significant (more explicitly, if there is at least one eigenvalue of 
A which is larger than unity). Therefore, our method should also have some limitations on its convergence and 
must be applied with care. 
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